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ABSTRACT 


Ahmed body is a standard configuration of road vehicles and most of the studies of automobile aerodynamics 
are performed based on it. In this paper, the plasma actuator was used as an active flow control method to 
control the flow around the rear part of the Ahmed body with the rear slant angle of 25°. Experiments were 
carried out in a wind tunnel at two different velocities of U=10m/s and U=20m/s using steady and unsteady 
excitations. The hot-wire anemometer was used to measure the vortex shedding frequency at the downstream 
of the body. Pressure distribution was measured using 52 sensors and total drag force was extracted with a 
load cell. Furthermore, smoke flow visualization was employed to investigate the flow pattern around the 
body. The results showed that the plasma actuator was more effective on the pressure distribution and total 
drag force at the velocity of U=10m/s. In fact, by applying steady and unsteady excitations there was 7.3% 
and 5% drag reduction; respectively. While at the velocity of U=20m/s; the actuator had no significant effect 
on pressure distribution and total drag. As a remarkable result, the plasma actuator, especially in the steady 
actuation, has demonstrated its effectiveness on dispersing the longitudinal vortices and suppressing the 
separated flow on the rear slant at low velocities. 


Keywords: Automotive; Aerodynamics; Wind tunnel; Drag reduction. 
NOMENCLATURE 


A effective frontal area Px static pressure of free flow 


AC altering current Pp power of plasma actuator 
Cp pressure coefficient pp pick to pick 

Cp drag coefficient Re Reynolds number 

D drag force TRMS True Root Mean Square 
H model height U free stream velocity 

In induced W model width 

L model length p density 

P pressure on the surface 


consumption and consequently the lower pollutant 
emissions (Hucho 1998). 


1. INTRODUCTION 


In recent years, by increasing the number of road 
transportation vehicles and consequently high fossil 
fuel consumption, the level of air pollution has 


Since the experimental study of flow control over 
different car models is expensive and time- 
consuming, normally a standard model is used as a 


raised to an alert situation in many metropolises 
(International Energy Agency 2007). In this regard, 
considering the international policies for reducing 
the consumption of fuels, the importance of the 
aerodynamics has been highlighted to minimize the 
resisting force of air flow. The lower amount of this 
force (known as drag force), causes lower fuel 


sample of all car models. Experimental tests are 
carried out on the standard model and results will be 
applied for designing real cars. A standard model 
known as Ahmed body (Ahmed ef al. 1984) which 
flow over it simulates three-dimensional flow 
pattern around an actual car due to its rear slant 
angle is popular in this regard. In some of the 


S. Shadmani et al. / JAFM, Vol. 11, No. 5, pp. 1267-1276, 2018. 


researches carried out on the rear part of the model, 
flow pattern was investigated by changing the rear 
slant angle (Tunay et al. 2014; Strachan et al. 2007; 
Vino et al. 2005; Tunay et al. 2016). The results of 
these studies showed that the flow pattern can be 
considered two-dimensional for a < 12.5°, and the 
drag force is independent of the angle variations 
(Strachan et al. 2007). For slant angles of 
12.5°<a<30°, the drag force increases with 
increasing the angle and the flow pattern is 3- 
dimensional (Fig. 1). The drag force reaches its 
maximum value at a=25° (Fig. 2) (Vino et al. 
2005). For a>30°, the flow is totally separated and 
the drag force is independent of the angle (Tunay ef 
al. 2016). Therefore, the flow regimes in the first 
and third a-ranges are similar and fall within low- 
drag zones. However, due to the three-dimensional 
flow, the second a-range (12.5°<a<30°) is known 
as high-drag zone (Fig. 2). 


In other works, active and passive flow controls 
were adopted in the tailing section of the model in 
an attempt to modify the flow in this region and 
reduce the drag force. Installation of various types 
of deflectors (Hanfeng et al. 2016; Fourrie et al. 
2011) is among the passive flow control methods 
which require no external energy (Hanfeng et al. 
2016). One of the drawbacks of these methods is 
that they have uncontrolled performance for various 
conditions (Kourta et al. 2012). On the other hand, 
active control methods such as mixed jets (Kourta 
et al. 2012), pulse jets (Joseph et al. 2012) and 
fluidic oscillators (Metka et al. 2015) have been 
investigated. Controllable nature of these methods 
allow them to be turned on/off on demand. Results 
of these methods in terms of reducing the drag force 
applied to Ahmed model are given in Table 1. 


Dielectric barrier discharge (DBD) plasma actuator 
is another modern tool for active flow control which 
has been used in aviation applications. Figure 3 
depicts the plasma actuator mechanism. It consists 
of two electrodes separated by a dielectric barrier, 
whereby air molecules above the insulated electrode 
are ionized by establishing a strong electric field. 
This ionized gas which is called plasma extends 
from the edge of the upper exposed electrode to the 
trailing edge of the lower insulated electrode. This 
adds up local momentum to the flow passing above 
the region as moving charged particles collide with 
other neutral particles of the gas surrounded the 
actuator (He et al. 2007). 


These actuators have many applications including 
separation control on leading edge of airfoils (Post 
and Corke 2004; Benard et al. 2008; Puoryoussefi 
et al. 2016), dynamic stall control of airfoils (Post 
and Corke 2006), flow control on bluff bodies (Do 
et al. 2007; Thomas et al. 2008; Thomas et al. 
2005; Rizzetta and Visbal 2008; Gregory et al. 
2008), flow control on boundary layer (Schatzman 
and Thomas 2008; Baughn et al. 2006; Font 2006), 
lifting applications (Corke et al. 2004), and flow 
control in turbomachines (Huang et al. 2006; 
Vanness et al. 2006). 

Boucinha et al. (2011) installed three groups of 
plasma actuators on the rear slant of an Ahmed 
body to investigate only the drag reduction on the 


model at velocity of 10m/s. They showed that the 
use of these actuators can reduce the drag force up 
to 8% by adopting the actuators on the upper part of 
the surface. Moreover, Khalighi et al. (2016) 
installed plasma actuators on the four edges of an 
Ahmed body model with 0° rear slant angle to 
analyze passing flow from the trailing part of the 
model. They showed the drag force decreases up to 
21.4% at 10m/s. 


In the present work, effects of plasma actuator 
system for flow control around an Ahmed body 1s 
experimentally studied in a wind tunnel using 
pressure distribution, drag force measurement and 
flow visualization. Experiments were performed at 
Reynolds numbers of Re=4.5x105 and Re=9x105 
(based on model length of 0.67m and free flow 
velocities of U=10m/s and U=20m/s, respectively), 
with the rear slant angle of 25° to have a dominant 
three-dimensional flow pattern and highest drag 
force (Fig. 2). Moreover, measurements of the 
natural frequency of vortex shedding, unsteady 
actuation, and flow control capabilities were 
compared for steady and unsteady actuations. In 
previous researches using all these methods and 
conditions has not been reported. It should be noted 
that the plasma actuator installation position may 
correspond to trailing edge of the roof, trunk edge, 
or below rear bumper of a vehicle. 
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Fig. 1. Flow pattern behind Ahmed body model 
at the rear slant angle of 25° (Ahmed et al. 1984). 
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Fig. 2. Plot of drag coefficient vs. slant angle 

(Ahmed et al. 1984). 


2. EXPERIMENTAL DETAILS 


Experiments were performed in an open loop low- 
speed wind tunnel at K. N. Toosi University (Iran) 
with a test section of 1x1.2m? and length of 3m 
(Fig. 4). Turbulence intensity of the flow was below 
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0.2% and non-uniformity of the flow along the test 
section was about +5%. 


2.1 Model Description 


The model was cons built of 6mm-thick Plexiglas at 
64% of the full scale original model. As mentioned 
before, in order to establish three-dimensional flow 
behind the model and maximize the drag force, the 
rear slant angle was 25°. Figure 5 shows the details 
of the model. The blockage ratio of the model in the 
test section is 3.8%. 


In order to eliminate the wall effects on the bottom 
of the wind tunnel (Hucho and Sovran 1993), the 
model was mounted on a rectangular surface at 
given spacing (Fig. 6). The distance between 
surface and tunnel floor was 16x10-2m, with the 
model connected to the surface via four cylindrical 
supports with the diameter of 2x10-2m and height 
of 32x10-2m. Moreover, leading edge of the 
surface was designed sharply to minimize the 
thickness of the developed boundary layer (Fig. 7). 
The thickness of the boundary layer formed on the 
surface at flow velocities of 10 to 50 m/s was 
about 1.5x10-2m (Joseph et al. 2012) so that the 
boundary layer on the surface has no effect on the 
model. 


2.2 Instruments and Measurement 


Techniques 


Pressure distribution was obtained using 52 
pressure sensors (pressure transducers) installed 
only on half of the slant and vertical surfaces of the 
model due to symmetry. The arrangement of the 
sensors is demonstrated in Fig. 8. Each sensor 
consists of a tube of outer and inner diameters of 
4mm and 2mm respectively. The pressure sensors 
were connected to a 28-channel pressure transducer 
apparatus (Honeywell-DCOO5NDC4). The results 
were monitored on a computer using an A/D 
(PCI6224) board, Lab View software, and pressure 
distribution software. Sampling interval was set to 5 
second and pressure coefficients were calculated 
using the Eq. (1): 


(1) 


Where P is pressure on the surface, Pæ is the 
static pressure of free stream, Ux is the velocity 
of free stream, and p is air density in the 
environment. The uncertainty of the 
measurement device is about +0.01. 


The drag force was determined using a load cell. In 
order to minimize friction force between the model 
and lower surface, the model was positioned on two 
rows of rails and wagons (Fig. 9a), and according to 
Fig. 9b, attaching a light aluminum rod of 0.5 m 
which was attached to load cell at +1 gr accuracy on 
the other end. 


In order to find natural frequency of vortex 
shedding behind the model, Hot-wire anemometer 
system was utilized (Fig. 10). The Hot-wire probe 
could be displaced along x and y directions with a 
traverse mechanism. 
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Table 1 Results of the research on the reduction 
of drag force. 
Drag 


Reynolds Reduction 


Number 


(%) 
Hanfeng et al. 5 
2016 8.7x10 9.3-10.7 


Fourrie et al. 2011 | 3.1-7.7x10° 


Studies 


Kourtaa et al. 2012 | _ 1.2x10° mes 
Joseph et al. 2012 | 1.4x10° | 68 | 
Metka et al.2015 | 1.4x10° 





Fig. 4. The used wind tunnel in this research. 
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Fig. 5. Three-dimensional view of Ahmed model 
and model dimensions (mm). 


m— 2W— We 2 





ma 2L -a L -= L m 





Fig. 7. View of the leading edge of the table 
which is designed sharply. 


2.3 Plasma Actuator System 


The input voltage to the plasma actuator was supplied 
by a high-voltage AC with the sinusoidal carrier 
wave and maximum output power of 1000W. Two 
digital multimeter sets were used to check out 
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actuation and carrier frequencies. The duty cycle and 
applied voltage were monitored by an oscilloscope. 
Furthermore, a digital multimeter (True RMS) was 
used to measure the mean electric current (Fig. 11). 


According to Fig. 12, the plasma actuator mounted 
such that, the exposed electrode was installed on the 
trailing edge of the upper surface, while the covered 
electrode was placed on the leading edge of the 
slant, in an edge-to-edge arrangement. The 
electrodes were made of copper, with the thickness 
of S5Omicrons. The length of the electrodes was 
equal to the model width, 25cm, with the widths of 
the exposed and covered electrodes of 5 and 15mm, 
respectively. In all arrangements, plasma actuator 
dielectric was composed of six layers of Kapton 
adhesive at the breakdown voltage of 7kV/mm and 
dielectric constant of 3.4. 


For a quiescent flow, installing a plasma actuator on 
a flat surface, a digital manometer (Testo 0560 
5126) along with a silicon micro-tube were applied 
to determine maximum induced velocity versus the 
applied voltage, Vpp  (pick-to-pick). The 
corresponding voltage was constant in all 
experiments. Other electrical parameters including 
carrier frequency (fac), was set to 10kHz, and 
actuation frequency, was determined based on the 
natural frequency of vortex shedding behind the 
model. Moreover, the actuator acted at steady and 
unsteady excitations, and duty cycle in unsteady 
actuation was 50%. 


3. RESULTS AND DISCUSSION 


3.1 Quiescent Flow Test 


In the first experiment, the maximum induced 
velocity in quiescent conditions was measured. Here 
we tried to reach the appropriate carrier voltage to 
reach the maximum induced velocity. The results are 
depicted in Fig. 13. It can be seen that the mean 
velocity occurs at 6kV. The induced velocity reached 
a maximum of 6m/s and then dropped down with 
further increase in voltage. This indicates that, with 
increasing the applied voltage to overcharged state, 
the plasma actuator becomes saturated, thereby 
changing the plasma regime from uniform mode 
(Fig. 14a) to stratified mode (Fig. 14b). 


In this situation, due to high conductivity, plasma is 
discharged at the higher density, reducing plasma 
density outside the region. This eventually leads to 
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Fig. 8. The position of pressure sensors. 


the reduction of induced velocity and this reduces 
actuator body force (Pouryoussefi et al. 2015; 
Dalvand et al. 2018). 






Fig. 9. a) Rail and wagon mechanism and; 
b) Load cell mechanism to measure drag force. 


Fig. 10. View of the position of Hot-wire probe 
behind the model. 








Fig. 12. The position of plasma electrodes. 
3.2 Hot-Wire Measurement 


A proper carrier frequency in unsteady excitation 1s 
estimated using vortex shedding frequency. In most 
of the time, this frequency is set to be equal with the 
shedding frequency (Pouryoussefi et al. 2015). For 
this purpose, 1-D Hot-wire probe was placed at the 
three different locations behind the model shown in 
Fig. 15. 


Figures 16-17 show spectral analysis of hot-wire 
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signal (effective instantaneous velocity) at (X°=0.5, 
Y"=0.6) position to demonstrate variations of 
velocity at the natural frequency of vortex shedding. 
At U=10m/s and U=20m/s, dominant frequencies of 
vortex shedding (peak of the curve) were about 
20Hz and 40Hz, respectively. Therefore, the 
dimensionless Strouhal number would be equal to 
0.37 at both velocities, which is in good agreement 
with other works reported in Table 2. Strouhal 
number is calculated based on Eq. (2): 





St = 





fxH (2) 
U Fig. 14. a) Uniform plasma distribution regime 


and; b) stratified plasma distribution regime 
Where H is the model height, f is natural frequency 


of back-model vortices, and U is free stream 
velocity. 


Considering the results obtained in this section and 
the section 3-1, the input electrical parameters 
which are used in all experiments and also the 
power consumption of the plasma actuator are 
presented in Table 3. The power of plasma actuator 5 
is calculated according to the research work of j 
Pouryoussefi et al. (2016) as Eq. (3): H 








P= VPP x I pis x COS (3) Fig. 15. Position of 1-D Hot-wire probe behind 
a the model. 
Where Jrus is ampere, Vpr is pick to pick voltage ‘sel f = 20Hz 


and ọ is about 6 ° (Pouryoussefi et al. 2016). 


6 


Amplitude 





Shell sind 


4 Frequency f (Hz) 





A Fig. 16. Plot of frequency and amplitude of the 
E vortexes behind the model at the velocity of 10m/s. 

2 

‘ E f = 40Hz 
0 ; \ 
0 2 4 6 8 Al nay 
Vpr (kV) | 

Fig. 13. Plot of induced velocity by plasma Frequency f (Hz) 
actuator, Um (induced flow) based on The Fig. 17. Plot of frequency and amplitude of the 

applied voltage, Vpp (pick-to-pick). vortexes behind the model at the velocity of 20m/s. 


Table 2 Comparison between the obtained Strouhal number in this research and other research works. 


Table 3 Input electrical parameters of plasma actuator. 


Veloci Carrier Voeltaee Power Excitation 
Oí aa Type Frequency y kv) Consumption Frequency 
fac (kHz) = Pp (W/cm) (Hz) 


aio lsg o e as i oa Tiot 





2 | 10 | Unsteady | 10 | 6 f 202 | 170 | 50 | 2 | 
3 | 20 | Steady { 10 {| 6 | 86 { 24) {too fs 
| 4 | 20 | Unsteady | 10 | 6 | 202 | 170 
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3.3 Pressure Distribution 


Pressure distribution difference between the front 
and rear of the model leads to pressure drag force. 
Krajnovic and Davidson (2005) showed that the 
contribution of pressure drag (Form drag) is about 
80% of total drag. This necessitates obtaining the 
pattern of pressure variations across two end 
surfaces of the model to investigate the effect of 
flow across the region. Figure 18 shows pressure 
coefficients contours on the slant surface behind the 
model at U=10m/s (Re=4.5x10°) and U=20m/s 
(Re=9x10°), which is good in agreement with the 
results reported by Keogh et al. (2016), Lienhart et 
al. (2003) and Joseph et al. (2012) for no actuator 
conditions (Fig. 19). 


According to Fig. 18, low-pressure regions of the 
flow extend along leading edge of the slant so that 
one can stipulate, separation bubble and flow jump 
over the leading edge of slant have occurred in such 
regions (see Fig. 1). Moving toward central areas, 
the low-pressure zones are gradually replaced by 
higher pressure, i.e. the pressure is recovered. At 
the lateral edges, low-pressure are developed on the 
surface, this indicates the isolation of C-pillar 
vortexes from the top corner of the surface along 
the surface length (see fig. 1). 


It can be observed in Fig. 20, by applying plasma 


Cp: O08 O76 07 066 06 055 05 04 04 14 D3 





Y (em) 


* (us) oz 


actuation at the velocity of U=10m/s, the volume of 
the low-pressure zones on the leading edge of the 
surface reduced considerably, such that the region 
was completely disappeared in steady state 
actuation. This indicates that by actuating the 
plasma, the developed vortex is degraded and the 
separated flow tends to attach to the surface. This 
will be observed in visualization tests as well. 
Moreover, reduction of low-pressure zones in the 
lateral edges of the surface 1s assign of the reduced 
effect of C-pillar longitudinal vortexes behind the 
model. The effect of steady actuation can be seen 
more significant than that of unsteady one. 


According to Fig. 21, at the velocity of U=20m/s, 
plasma actuation imposed no significant effect on 
the reduction of low-pressure regions developed by 
the vortices behind the model, such that nearly the 
same pressure distribution is obtained for the steady 
and unsteady actuations. 


Figure 22 shows the comparison between pressure 
distribution of no-actuation, steady actuation, and 
unsteady actuation conditions on the centerline of 
the model at velocities of 10m/s and 20m/s. At the 
velocity of 10m/s, plasma actuation has resulted in 
considerable rising pressure coefficients, 
particularly across the regions at the slant while the 
effect of plasma actuator is insignificant at the 
velocity of 20m/s. 





Fig. 18. Contours of pressure distribution on the rear slant of an Ahmed body in plasma off mode at 
velocities of a) 10m/s and; b) 20m/s. 
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Fig. 19. Contours of pressure distribution on the rear slant of the Ahmed model in velocity of 20m/s 
and no actuation condition in the studies of a) Keogh et al. (2016); 
b) Lienhart et al. (2003); c) Joseph et al. (2012) and; d) Present study. 
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Fig. 20. Contours of pressure coefficient distribution on the rear slant of Ahmed body at the velocity of 
10m/s in a) Plasma off; b) Steady actuation and; c) Unsteady actuation. 
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Fig. 21. Contours of pressure coefficient distribution on the rear slant of Ahmed body at the velocity of 
20m/s in a) Plasma off; b) Steady actuation and; c) Unsteady actuation. 
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Fig. 22. Plot of pressure coefficient distribution on the centerline 
of the model at a) U=10m/s and; b) U=20m/s. 


3.4 Drag Measurement 


In the present research, total drag force was 
measured using a load cell. For no plasma 
actuation, the drag force was found to be 0.8N at 
Re=4.5x105 (U=10m/s) and 3.22N at Re=9x105 
(U=20m/s). Since drag force 1s directly proportional 
to the velocity squared, it can be observed that 
doubling the velocity, leads to four times higher 
drag force. According to the drag force coefficient 


(Eq. (4)): 
Fp 


Cy =———— 
0.5x px AxU? 


(4) 


Where U is free flow (stream) velocity, A is the 
effective frontal area, and p is air density. 


The drag coefficient of the model was found to be 
0.29, which is in good agreement with the results 
reported by Ahmed et al. 1984 (see Fig. 2). 


Based on Fig. 23, applying plasma actuation at the 
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velocity of 10m/s, the values of drag force reduced 
about 7.3% and 4.88% in steady and unsteady 
actuations, respectively. While, at the velocity of 
20m/s, the reduction of drag force in steady and 
unsteady actuations reached 2.13% and 1.5%, 
respectively. The results indicate that, at lower 
velocities, particularly in steady actuation, the 
plasma actuator can reduce overall drag force more 
effective. Considering the fact that, the dominant 
drag force in the Ahmed model is pressure drag 
(Krajnovic and Davidson 2005b), the results 
obtained from drag force measurement using load 
cell shows good agreement with the pressure 
measurement results. 


3.5 Flow Visualization 


Since the plasma actuator was more effective at 
10m/s, for flow visualization study we only 
conducted our experiment at this velocity. 
According to Fig. 24 a, in plasma off mode, the 
airflow separates from the surface as it passes over 
the leading edge of the rear slant, developing a 
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vortex on the slant surface. Moreover, within the 
flow passing through the lower part of the model, a 
larger vortex 1s formed behind the vertical surface. 
Figures 24b and 24c show that the plasma actuator 
suppresses the flow passing over the surface edge 
by actuating shear layer, and the separated flow 
from the slant surface gets closer to the surface. As 
a result, the developed vortex on the slant is 
degraded and the volume of the wake formed 
behind the vertical surface is reduced. This lead 
reduction in pressure drag (Form drag) (Section 
3.4). 


8 


6 


Drag reduction (%0) 
hes, 






8t.l0m/s 8 UnstlOm/s 88t20m/s m Unst20m’s 
Fig. 23. Plot of drag reduction 
(St: Steady actuation, Unst: Unsteady actuation). 
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+ ofa > 
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Fig. 24. Pattern of the flow passing through the 
region behind the model at U=10m/s in: 
a) Plasma off; b) Steady actuation and; 
c) Unsteady actuation. 


4. CONCLUSIONS 


In the present research, the plasma actuator was 
used to control flow over an Ahmed body model 
with a rear slant angle of 25°. Plasma actuator 
electrodes were installed on the leading edge of the 
slant behind the model, and actuations were 
performed in steady and unsteady excitations at 
different Reynolds numbers with flow velocities of 
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U=10m/s and U=20m/s. 


In quiescent flow test, using the plasma actuator on 
a flat surface, the appropriate applied voltage was 
determined based on maximum induced velocity. 
Also, actuation frequency of the plasma actuator 
was obtained by measuring natural frequency of the 
vortex shedding behind the model at U=10m/s and 
U=20m/s velocities. 


Visualization results indicate that at U=10m/s the 
plasma actuator is capable of suppressing the 
separated flow from the leading edge of the rear 
slant and attaching the flow back to the surface by 
actuating the shear layer. For this case, the volume 
of longitudinal vortices behind the model is reduced 
while increasing pressure coefficient on the slant at 
the same time. By reduction of the pressure 
difference between front and rear part of the model, 
pressure drag force and hence total drag force 
decreases, so that the reduction of drag force in 
steady and unsteady excitations was 7.3% and 
4.88%, respectively. At the velocity of U=20m/s, 
the results indicated that the plasma actuator is not 
effective on pressure enhancement on the rear part 
of the model and the reduction of drag force in both 
steady and unsteady excitations is about 2%. 
Therefore, application of the plasma actuator at low 
velocities, particularly in steady excitation, can 
serve as an effective tool for flow control and drag 
reduction. 


For further investigation of plasma actuator 
effectiveness, it is suggested to use several plasma 
actuators in different positions of the model 
simultaneously and evaluate the relationship 
between the plasma actuator energy consumption 
and the amount of fuel reduction caused by the drag 
force reduction evaluated. 
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